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Abstract 

The introduction of a metric onto the space of parameters in models in Statistical 
Mechanics and beyond gives an alternative perspective on their phase structure. 
In such a geometrization, the scalar curvature, 1Z, plays a central role. A non- 
interacting model has a flat geometry (1Z = 0), while 1Z diverges at the critical 
point of an interacting one. Here, the information geometry is studied for a number 
of solvable statistical-mechanical models. 
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1 Introduction 



As the application of statistical-physics techniques becomes more widespread 
and acceptable outside the traditional physics community, it is clear that the 
phenomena of phase transitions have important roles to play there. Indeed, 
phase transitions are common to a very wide range of disciplines, from physics 
to biology, economics and even to sociology. As statistical physics finds perti- 
nence in these areas, so too can the field draw on concepts outside the confines 
of pure physics. 

In all of these disciplines, models are characterised by certain sets of param- 
eters. The idea of endowing the space of such parameters with a metric and 
geometrical structure has been borrowed from parametric statistics [1]. Given 
a probability distribution p(x\9), and a sample X\, . . . , x n , the objective is to 
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estimate the parameter 9. This may be done by maximizing the so-called like- 
lihood function, L(9) = Yl7=iP( x i\^)j 01 its logarithm (called the log-likelihood 
function), 



lnL(0) =£>pW0) 



(1) 



i=l 



The gradient of this quantity is the score function: 
d\nL(9) 



U{9) 



d9 
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and the expectation of this random variable is zero (E[[/(0)] = 0). Its variance 



is 



Vai[U(6)] = E 



-d 2 \nL(6) 
d9 2 



(3) 



This quantity is called the expected or Fisher information. Taylor expanding 
the log-likelihood function, one arrives at 



kiL(0 + 66)-lnL(0) = 59 
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dO 



+ 



(59) 2 d 2 \nL(9) 



d9 2 
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The first term on the right-hand side is zero at the true #-value. Therefore, 
the closeness of two probability distributions characterised by 9 and 9 + 59, is 
given by the second term, or the Fisher information. 

For higher- dimensional distributions (which may be continuous), where, in- 
stead of 9, one has a set of parameters, 9i,9 2 , ■ ■ ., the Fisher information is 
defined as 



GiM = -E 



d 2 \np(x\9) 
89 i d9 j 



p(x\9) 



d 2 \np(x\9) 
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dx . 
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C.R. Rao suggested this is a metric. It is, in fact, the only suitable metric in 
parametric statistics and is called the Fisher- Rao metric [1]. 

In generic statistical-physics models, we have two parameters, f3, which we 
may think of as the inverse temperature, and h, the external field. In this case 
the Fisher-Rao metric is simply given by 

Gij = didjf , (6) 
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where / is the reduced free energy per site and = (d/d/3, d/dh). 
For such a metric the scalar curvature may be calculated as 



K = 



2G 2 



djf d p d h f dlf 

a?/ dp h fd p dif 



(7) 



where G = det(Gjj) is the determinant of the metric itself. The scalar curva- 
ture plays a central role in any attempt to look at phase transitions from a 
geometrical perspective. Indeed, 1Z measures the complexity of the system. A 
flat metric implies that the system is not interacting. Conversely, and for all 
the models that have been considered so far, the curvature diverges at (and 
only at) a phase transition point for physical ranges of the parameter values. 

Using standard scaling assumptions, we can anticipate the behaviour of 1Z 
near a second-order critical point. With t — 1 — (3/(3 c , 

f((3, h) = A" 1 / (t\ at , h\ a ») = tk^ (ht'^j , (8) 



where 



1 ^ (Q\ 

at = Vd ' ah = Vd ' (9) 



are the scaling dimensions for the energy and spin operators and d is the 
spacial dimensionality. One finds for the scalar curvature, 
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(10) 



and 



G ~ t a t 



+2 



yielding 



n ~ r ~ \P - A 



io-2 



:i2) 
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where hyperscaling {yd = 2 — a) is assumed and £ is the correlation length. 

The systems hitherto analysed from the information-geometry perspective in- 
clude the Ising model in one dimension [2], the Bethe lattice Ising and mean- 
field models [3]. The main results hitherto established are that 1Z is positive 
definite and diverges (as £ d ) only at the critical point. 

In an effort to see how generic these features are, and to discover new ones, we 
analyse the Potts model in one dimension, the Ising model in two dimensions 
coupled to quantum gravity and the spherical model in three dimensions. 



2 Information Geometry in Specific Models 



The Potts Model in One Dimension: The one-dimensional g-state Potts 
model, like its Ising counterpart (which corresponds to q — 2) is exactly solv- 
able. Although it has no true phase transition, thermodynamic quantities di- 
verge at zero temperature. In the Ising case, explicit calculations showed that 

_ cosh h , . 

Rising = 1+ , • (13) 

/sinh 2 h + exp (-4/3) 



The correlation length is known to behave as £ ~ — 1/ln (tanh/3), and there- 
fore 7^isi ng ~ £ as (3 — > oo. 

For plotting purposes, it is convenient to introduce 

y = exp (/?) and z = exp (h) . (14) 



The curvature is plotted in Fig. 1. It is clear that 7^i S i ng is positive definite and 
exhibits a h — > —h (or 2 ^ 1/z) symmetry. Furthermore, given a particular 
value for y, TZj S i ng is maximum along the zero-field ridge z — 1. 

The one- dimensional g-state Potts model is also exactly solvable, so it is sen- 
sible to exploit it in an attempt to decide which of the previous features are 
generic. An explicit calculation [4] shows that 

T^Potts = A(q, y,z) + ^ ~ y ~ £ , (15) 



for h = and y — > 00. Here, i](q, y, z) is the Potts analogue of the Ising term, 
sinh 2 h + exp (—4/3), and A and B are smooth functions of y and z. So the 
expected scaling (12) holds with d — 1 or a — 1. 
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Figure 1 also shows that, in the Potts case, 7?potts is no longer positive definite. 
Furthermore, the z —> l/z symmetry of the Ising model is no longer present. 
For a given y, 7Zp otts no longer peaks along z — 1. Finally, an analysis of 
the Lee- Yang (complex h) zeroes of the model reveals that the curvature also 
diverges as the locus of zeroes, zo(y), is approached. In fact, 

K~(z-zo(y))-* , (16) 



so that this divergence is characterised by an edge exponent a = 1/2. 




The Three-Dimensional Spherical Model (and the Ising Model on 
Planar Random Graphs): The spherical model is given by 



J spherical 



J ds\ . . . ds L d5 ^2 s 1 ~~ exp 



L (ij) 



, (17) 



and can be solved by exponentiating the constraint and using steepest descent. 
Allowing the lattice extent, L, to diverge reveals no transition for d — 1,2. 
The transition for d = 3 has a = —1, (3 = 1/2 and 7 = 2. This is the 
same set of critical exponents as for the Ising model on planar random graphs 
(matter coupled to 2D gravity). This is remarkable, because there are no 
obvious similarities between the two models. 

Explicit calculations (in both models) give [5,6] 

n~\f3-(3 c \-* , (18) 
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which does not accord with the prediction of (12). That prediction is 1Z ~ \/3 — 
(3 c \ a ~ 2 , with a — — 1. The source of the discrepancy can be traced back to the 
top left term in the determinant in (10), which is t~ a . In both current models, 
a is, in fact, negative and this term vanishes as criticality is approached. It is 
replaced by a constant term coming from the regular part of the free energy. 
Both models then yield the same result (18). 

3 Conclusions 

In statistical physics, and in related fields - from the bio-sciences to economics 
- phase transitions play a central role. Thus new insights into the character- 
isation of critical phenomena are of paramount importance. Here, geometric 
ideas from the field of parametric statistics are "borrowed" and explored. It 
is found that the curvature associated with the Fisher-Rao metric is, indeed, 
a useful quantity in the characterisation of phase transitions. Some features 
found in the Ising model are found not to be generic. In particular, TZ can be 
negative and there is no symmetry nor ridge along h = 0. More surprisingly, 
the naively expected scaling behaviour (12) fails in both the three-dimensional 
spherical model and in the Ising model on planar random graphs. Reasons for 
this are given. Once again, it is curious that all critical exponents for the latter 
two models coincide, although there is no obvious physical relation between 
them. 
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